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We examine two-spin (2S) and cyclic four-spin (4S) magnetic interactions 
in insulating ladder and two-dimensional (2D) cuprates. By a comparison 
^ ■ of eigenstates between d-p and Heisenberg models, we evaluate magnitudes 

O I of these interactions. We find that the 4S interaction is ^10 % of near- 

est neighboring 2S interaction, and a diagonal 2S interaction is considerably 
small. The 4S interaction for a ladder cuprate is larger than that for 2D 
^ , one, and plays an important role in the low- energy excitation. The Heisen- 

'nI" I berg ladder with the obtained 2S and 4S interactions reproduces very well the 

experimental result of the temperature dependence of the magnetic suscepti- 
\^ . bility. 

O' PA CS numbers: 71.70. Gm, 75.30.Et 

C^ ■ A variety of insulating cuprates, including the parent compounds of 

high- Tc superconductors, afford us an opportunity to study magnetic prop- 
,.^ , erties of low-dimensional systems. They have unique Cu-0 structures formed 

j;H ' by CUO4 tetragons, that is, two dimensional (2D) Cu02 planes, two-leg lad- 

O ■ ders and one dimensional (ID) chains. In insulating cuprates, magnetic 

interactions J between Cu S=l/2 spins become important key parameters 
for understanding of the electronic states. In 180° Cu-O-Cu coupling, which 
^ ■ is realized when CUO4 tetragons share their corners, the superexchange pro- 

cess via O ions gives a dominant contribution to J. 

Recent magnetic measurements demonstrate remarkable dependence of 
J between nearest neighboring (NN) Cu spins on the dimensionality of Cu- 
O structures (i) J in ID cuprates is larger than that in 2D ones, and (ii) 
for ladder cuprates, J along leg direction is larger than that along rung 
one. These can not be understood only by considering the bond length 



dependence of J. In our previous studyJJ we have shown that the hopping 
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matrix elements between Cu3d and 02p orbitals (Tpd) and between 02p 
ones (Tpp) are considerably influenced by Madelung potentials around Cu-0 
and 0-0 bonds. As a result, J depends on the crystal structure through 
the modification of Tp^ and Tpp due to the Madelung potential. By taking 
into account these effects, we have clarified the microscopic origin of the 
characteristics of (i) and (ii). 

Very recently, an inelastic neutron experiment for the two-leg ladder 
cuprate La6Ca8Cu2404i has been performed in the wide range of the excita- 
tion energy.D By comparing the observed dispersion of spin excitation with 
the theoretical one, it has been shown that J along leg direction is larger that 
along rung one, and in addition a four-spin (4S) interaction (~0.005 eV) is 
necessary to explain the experiment. Here, the 4S interaction is given by the 
cyclic permutation on four Cu spins forming a plaquette (see Fig. l(a-2)). 
This result indicates importance of not only usual two-spin (2S) interac- 
tions but also 4S interaction among Cu spins in a plaquette. Therefore, it is 
necessary to establish a proper magnetic description for undoped cuprates. 
This will also yield an appropriate understanding of magnetic properties in 
high- Tc cuprates. 

So far, some theoretical studies on the 4S interaction have been done, 
but they are only for 2D cupratesBo In this paper, focusing on insulating 
ladder and 2D cuprates, we investigate possible magnetic interactions in- 
cluding 4S one in the systems. By associating a d-p model with an effective 
Heisenberg model, we evaluate magnitudes of these magnetic interactions. 
The effect of 4S interaction on the magnetic properties is examined through 
the magnetic susceptibility, and is shown to lead to a dramatic change in 
the low-energy excitation for ladder cuprates. 

As a starting model for a ladder cuprate, we consider a d-p model 
as shown in Fig. l(a-l), where three and two Cu3d^2_y2 orbitals are ar- 
ranged along leg (x) and rung (y) directions, respectively, and 02po- orbitals 
(a = x,y) form the CUO4 tetragons (CugOir cluster). In the model, open 
boundary condition is imposed. The Hamiltonian to describe the model is 
expressed by Tpd and Tpp, the energy-level separation between Cu3d and 02p 
orbitals (A) and Coulomb interactions such as on-site Coulomb interactions 
at Cu site (Ud) and O one (Up). 

An effective model to describe the low-energy excitation of the d-p model 
can be given as a S=l/2 Heisenberg model with 2S interactions along leg 
(Jieg), rung (Jrung) and diagonal (Jdiag) and cyclic 4S interaction (Jcyc) as 
shown in Fig. l(a-2). The spin Hamiltonian is written by 



H — / ^ ('^IcgJi ■ ^i+x + '-'rung^i " ^i+y + •^diag^i " ^i+x+yj 
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+ Jcyc E (^Jkl + ^7kl)> (1) 

plaquette 

where Si is a spin operator at i site, and x (y) is a unit vector along x (y) 
directions. Here, Jcyc is defined as a coefficient of 4S cyclic permutation 
operators i-\jki and -Pi^i^i, which can be rewritten by using the 2S and 4S 
interactions as 

^jki + ^7ki = 4[(Si • Sj)(Sk • Si) + (Si • Si)(Sj . Sk) - (Si • Sk)(Sj • Si)] 
+ [(Si • Sj) + (Sj • Sk) + (Sk • Si) + (Si • Si)] + [(Si • Sk) + (Sj • Si)] + \. (2) 

We calculate the eigenstates and eigenvalues for the d-p model by nu- 
merically diagonalization method. For the ladder cuprate, SrCu203, Tp^^s 
are 1.08, 1.22 and 1.25 eV for leg, rung and inter ladder directions, respec- 
tively. A=2.6 eV, f7rf=8.5 eV, Up=4.1 eV and the direct exchange interac- 
tion between Cu3d and 02p orbitals (Kp^) is 0.05 eV. TppS are 0.52 eV and 
0.42 eV for solid and dotted bonds in Fig. 1(a), respectively.EJ Note that the 
magnitudes of two kind of Tpp's are different because of a difference of local 
environment, El i.e. the former has only one neighbor Cu sites, while the latter 
does two neighbor Cu sites, where one of them belongs to the neighboring 
ladder. 

By associating the obtained eigenstates and eigenvalues of the d-p model 
with those of the Heisenberg model of Eq. (1), we determine Jieg, Jmng, Jdiag 
and Jcyc. The five energy differences obtained by six lowest eigenvalues of 
the d-p model are used to fix these interactions. In addition to Cu^On 
cluster, we also used a quadratic CU4O12 cluster to determine the magnetic 
interactions for the 2D cuprate, La2Cu04 (see Ref. 1 for the parameters used 
in the calculation), as well as ladder cuprates. 

Table 1 

The magnetic interactions (energy unit in eV) for CU4O12 cluster (A) and 
CugOiy cluster (B) obtained by the fit to Heisenberg models. The numbers 
in parentheses represent the deviation in the last significant digit. 



2D cuprate (La2Cu04) 

t'NN «^diag "^cyc 



0.146(1) 0.000(0) 0.0108(8) 



Ladder cuprate (SrCu203) 

<^leg t^rung "^diag ^^cyc 



(A) 
(B) 



0.165(5) 0.15(0) 0.001(1) 0.015(0) 
0.195(5) 0.15(2) 0.003(2) 0.018(2) 



The results are summarized in Table 1. The results for a 2D cuprate 
are comparable with those by Schmidt et al..u Jcyc is ~7 % of the 2S NN 
interaction JnNj and Jdiag is zero. Jnn is ~ 0.15 eV, similar to our previous 
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study.EJ As for a ladder cuprate, we find that Jieg > Jrung in both CU4O 
and Cu60i7 clusters, which is again consistent with our previous study. 
The anisotropy is caused by the enhancement of interladder Tpp denoted by 
solid line in Fig. l(a-l). In addition, the enhancement enlarges the value 
of Jcyc for a ladder cuprate as compared with that for a 2D cuprate, and 
its magnitude is ~10 % of Jicg. Jdiag is almost zero. The energy levels of 
the Heisenberg model are shown in Fig. 1 (b) together with those for the 
d-p model, where all the eigenvalues Ei (i=0~6) of the Heisenberg model 
is shifted so that the energy of ground state Eq is equal to that of the d-p 
model. We can see a good agreement of the energy levels and eigenstates 
between both models. In order to check validity of the obtained interactions, 
we have also performed calculations for a d-p model with periodic bound- 
ary condition (CueOis cluster) and the corresponding periodic Heisenberg 
model. We confirmed that the Heisenberg model with the interactions in Ta- 
ble 1 reproduces very well the distribution of the eigenstates for the periodic 
d-p model. 

By using the magnetic interactions obtained, we examine the temper- 
ature dependence of the magnetic susceptibility x(^)- The exact diagonal- 
ization method is used to get the x(2^)'s for Heisenberg ladders with 2x4, 
2x6 and 2x8 sites. The results are shown in Fig. 3 for the parameters 
of Jieg=0.19 eV, Ji.ung=0.15 eV, Jdiag=0.003 eV and Jcyc=0.0175 eV and 
of the same values but Jcyc=0 eV. In the inset, x(^)'s at wider range of 
temperature are shown. 

The effects of Jcyc on x(^) are (i) to enhance the magnitude of x(^) at 
low temperatures and (ii) to shift a hump structure (~2000 K) in the lower 
temperature region. This suggests the decrease of an energy gap between 
the singlet and triplet states (spin gap). The decrease of the spin gap due 
to Jcyc is also observed in the dynamical spin correlation function.B The 
reason for the decrease is that Jcyc disturbs the singlet states on rungs as 
magnetic frustration, consequently making Jrung effectively decreased. We 
find that by considering Jcyc, x(^) is in agreement with the experimental 
xiT) for SrCu203.l3 Thus, it is interesting that Jcyc strongly influences the 
low-energy excitation above the spin gap. This is in contrast to 2D cuprates, 
where Jcyc only modifies the values of 2S interactions. Further study on the 
effect of Jcyc on other physical quantities is necessary for a quantitative 
investigation between theory and experiment. 

This work was supported by CREST and NEDO. The parts of the 
numerical calculation were performed in the Supercomputer Center in ISSP, 
Univ. of Tokyo, and the supercomputing facilities in IMR, Tohoku Univ. 
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Fig. 1. (a-1) The d-p model (CugOiy cluster) simulating a ladder cuprate. 
Unconnected orbitals denote Cu sites on the neighboring ladders. There are 
two inequivalent bonds for Tpp, which are distinguished by the solid and dot- 
ted lines, (a-2) The corresponding Heisenberg model, (b) The distribution 
of six lowest eigenstates for the d-p model (a-1) and the Heisenberg model 
(a-2). Indices in figure denote the irreducible representation and the total 
spin of the eigenstates. 
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Fig. 2. The calculated x(2^)'s for Heisenberg ladders. Experimental results 
are shown by filled circles.^ 



